Microtubule length dependence of motor traffic in cells 
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In living cells, motor proteins, such as kinesin and dynein can move processively along microtubule 
(MT), and also detach from or attach to MT stochastically. Experiments have found that, the traffic 
of motor might be jammed, and various theoretical models have been designed to understand this 
traffic jam phenomenon. But previous studies mainly focus on motor attachment/detachment rate 
dependent properties. Recent experiment of Leduc et al. found that the traffic jam formation of 
motor protein kinesin depends also on the length of MT [Proc. Natl. Acad. Sci. U.S.A. 109, 6100- 
6105 (2012)]. In this study, the MT length dependent properties of motor traffic will be analyzed. 
We found that MT length has one critical value N c , traffic jam occurs only when MT length N > N c . 
The jammed length of MT increases with total MT length, while the non-jammed MT length might 
not change monotonically with the total MT length. The critical value N c increases with motor 
detachment rate from MT, but decreases with motor attachment rate to MT. 



I. INTRODUCTION 

Motor proteins are essential for biophysical function- 
ing of living cells [H-Ql. Processive motors, such as ki- 
nesin and dynein, can transport organelles and vesicles 
along microtubule (MT) unidirectionally fsUlol]. Exper- 
iments found that the traffic of motors along MT might 
be jammed ll|, and actually, this jammed phenomenon 



existence and properties of domain wall (DW), which is 
defined as the interface of high motor density and low 
motor density 0, [l6 , 



21 



26|. 



One surprising finding in the recent experiment of 
Leduc et al. is that, different from macroscale vehicle 
traffic, the formation of traffic jam of motor protein ki- 
nesin (Kip3) depends not only on rates a, /3,uj a ,u)d but 



also on the length of MT 27|, l28[. For given rate con- 



has been predicted by various theoretical models |12h17| 



This traffic problem can be described by totally asym- 
metric simple exclusion process (TASEP), in which the 
track (MT) of motors is regarded as a one-dimensional 
lattice of size N. Each lattice site can be either empty 
or occupied by one motor. In the bulk, motors can hop 
from site i to site i + 1 with constant rate v (for sim- 
plicity, v = 1 is usually used), provided the target site is 
empty. At site 1, motors can enter the lattice from cell 
environment with density a , provided the site is empty. 
At site N, motors can leave the lattice into environment 
of density (3 with rate 1 — j3. Meanwhile, motors can at- 
tach to bulk sites of MT with rate oj a , and detach from 
bulk site of MT with rate u d . The TASEP has been 
studied extensively However, previous 

studies about TASEP mainly focus on rates a,f3,uj a ,ujd 
dependent properties in long MT length limit, such as the 



stants a, fi,u> a ,u)d, there exists one critical value N c of 
MT length, DW occurs only if the MT length N > N c . 
If there is no DW for N — > 00, there will be no DW for 
any finite value N. But, even if DW occurs for N — > 00, 
there will be no DW for MT length N < N c . In cells, MT 
length N is always finite, so it is biophysical necessary 
to study MT length N dependent properties of motor 
traffic. Since the rates a, /?, u) a , om dependent properties 



have been extensively analyzed [lj, ll6|, 1 1 81 - 4201 . |25I | , only 
the MT length N dependent properties of motor traffic, 
especially the rates a, (3, ui a , LUd dependent critical MT 
length N c , will be discussed in this study. 



II. METHOD 

Let pi be the probability that MT site i is occupied by 
one motor, then pi is governed by the following equations 
16| 
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dpi/dt = vpi-i(l— pi) - vpi(l- p i+ i)+u) a (l- pi) -UdPi, 

(1) 
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for 2 < i < N — 1, and at the boundaries i = l,N, 

dpi/dt = a(l-pi)-vpi(l-p2), ^ 
dp N /dt = vp N -i(l - p N ) - fip N . 

For kinesin Kip3 studied in [2^ . i = 1 is the minus-end 
of MT, and i = N is the plus-end of MT. For N large 
enough, Eq. ([I]) can be approximated as the following 
differential equation 

ep"{x) + v(2p(x) - l)p'(x) = (fi„ + n d )p(x) - fi„, (3) 

with boundary conditions p(0) = a/u, p(l) = 1 — /3/v. 
Where Q a ,d — lim]Vu ai ij, and e — >• is a small parameter 
[T4I ]. For the special case v = 1, the properties of mo- 
tor traffic, especially the existence of DW and boundary 
layer (BL) of motor density p{x) , have been detailed dis- 
cussed, see Tab. Uor [ljj. Let p; Q and be solutions of 
(2/o(x) — l)p'(x) — (fl a + Qd)p(x) — £l a but with boundary 
conditions pi a (0) = a and p r /s(l) = /3 respectively. The 
results listed in Tab. U can be roughly summarized as 
follows, (I) if < a < 0.5 and there exists Xo such that 
Pia(%o) + prpixo) = 1, then (a) DW occurs at x = xq 
if < xq < 1, (b) BL occurs at left boundary x = if 
xq < 0, and (c) BL occurs at right boundary x = 1 if 
xo > 1. Meanwhile, BL always occurs at right boundary 
x = lif0.5</3<l. (II) Otherwise, BL always oc- 
curs at left boundary x = 0, and will also occur at right 
boundary x = 1 if 0.5 < j3 < 1. 

However, as mentioned before, recent experiment of 
Leduc et al. found that motor traffic depends also on 
MT length N, DW occurs only when MT length N large 
enough [27]]. Therefore, the mean field results listed in 
Tab. U might not be reasonable. The MT length N de- 
pendent properties of motor traffic should be discussed. 
The results listed in Tab. U indicate that, for N — > 00, 
one necessary condition of DW formation is < a < 0.5. 
So in this study, we will mainly focus on this special case, 
and discuss the critical MT length N c . 

III. RESULTS 



the formation of DW and the ratio of jammed MT length 
dependent on MT length N. For 1/(K + 1) < J3 < 1, the 
ratio of jammed MT length increases with MT length N 
[see Figs. [[Ja, c)], but for < fi < 1/(K + 1) this ratio 
does not change monotonically with N [see Fig. [TJb)]. 
The results plotted in Fig. [TJd) show that the formation 
of BL also depends on MT length N. Let N L and N H 
be the length of MT with low and high motor densities 
respectively (in calculations, low density means motor 
density p < 0.5, and high density means p > 0.5). The 
results plotted in Fig. [5](a-d) show that, for the cases (2) 
(3), both the high density length Nh and its ratio Nh/N 
increases with MT length N, the low density length Nl 
first increases and then decreases with MT length N, but 
its ratio Nl/N decreases monotonically with MT length 
N. However, for the case (1), neither Nh/N nor N^/N 
changes monotonically with MT length N, though Njj 
and Nl increase with N for almost everywhere, see Fig. 
121(e) (f). 

One meaningful question about the MT length depen- 
dent traffic of motors is that when DW will occur, pro- 
vided DW occurs in large N limit. Let iV c be the crit- 
ical value of MT length, i.e. DW occures iff MT length 
N > N c . The curves plotted in Fig. [3] show that, N c 
decreases with motor attachment rate a to the first site 
of MT, but increases with detachment rate j3 from the 
last site of MT. Meanwhile, N c increases with attach- 
ment rate Lo a and deceases with detachment rate ujd- This 
means that, motors will be more likely to be jammed for 
high motor attachment rates (a and u>d) and low motor 
detachment rates (/3 and u>d)- 

Finally, the curves plotted in Fig. @]show that both the 
non-jammed MT length Nl and its ratio Nl/N increase 
with motor velocity v, while the jammed MT length Nh 
and its ratio Nh/N decrease with v. So slow motors are 
more likely to be jammed. One can easily deduce that the 
critical length N c of MT increases with motor velocity v, 
which means that long MT is needed to block the traffic 
of fast motors. 



The results of long MT limit listed in Tab. U indicate 
that, there are three typical cases in which DW might 
occur, (1) < a < 0.5, < /3 < l/(K + f), where K = 
uj a /uj d , (2) < a < 0.5, l/(K + I ) < f3 < 0.5, and (3) 
< a < 0.5 0.5 < f3 < 1 . The motor densities plotted in 
Figs. [IJa, b, c) show that, for any of the three cases, both 



IV. CONCLUSIONS 

Since the motor traffic along microtubule (MT) de- 
pends not only on motor attachment/detachment rates 
{lti| , but also on MT length 2jj| , the previous results ob- 
tained by assuming MT length N tends to infinity might 
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be misleading. Therefore, in this study, the MT length 
dependent properties of motor traffic are studied. We 
found that, for given attachment and detachment rates 
the MT length has one critical value N c , domain wall 
occurs only when MT length N > N c . This critical 
length N c increases with detachment rate but decreases 
with attachment rate. The jammed MT length Nu in- 
creases with MT length N, but depending on the mo- 
tor detachment rate f3 from the last MT site, the non- 
jammed MT length Nl might not change monotonically 
with N. Meanwhile, for not too small motor detachment 
rate, (3 > 1/(K + 1), the ratio of jammed MT length 
increases monotonically with MT length N. Calcula- 
tions also show that, both the non-jammed MT length 
Nl and its ratio N^/N increase with motor velocity v. 
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TABLE I: Rates a, p, cj a , Ud dependent properties of motor density p along MT for MT length N — > oo. Where a is the motor 
attachment rate to the first site of MT, /3 is motor detachment rate from the last site of MT, u) a ,oJd are the motor attachment and 
detachment rates to and from MT bulk sites, K = uj a /u)d- pi a , p r p, and p r are solutions of (2p(x) — l)p'(x) = (Q a +Qd)p(x) — fl a 
but with boundary conditions pi a (0) = a, prpiX) = P, and p r (l) = 0.5,p r (0) > 0.5, respectively. For further properties of 
DW, see The results listed here are only right for K = w a /wd > 1 and v = 1, results for K < 1 can be obtained by the 
motor-empty symmetry, and results for v ^ 1 can be obtained by replacing a, /3, ui a , Ud with a/v, fi/v, UJ a /v, Ud/v respectively. 
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FIG. 1: Four typical examples of motor density along MT with different length. For the sake of comparison, the MT length 
is normalized, (a) a = 0.24, P = 0.48, u a = 0.001, w d = 0.0002. Motor traffic is not jammed for short MT, but for MT length 
N > 120, motor is jammed from the right boundary. The ratio of jammed MT length Nh/N increases with MT length N. (b) 
a = 0.01, /? = 0.1, u a = 0.001, ujd = 0.0002. Similarly, motor traffic is not jammed for short MT, but as MT length N increases, 
motor is jammed from the right boundary. However, the ratio of jammed MT length Nh/N does not change monotonically 
with MT length N. (c) a = 0.25, /3 = 0.7, uj a = 0.001, w d = 0.0002. Similar as in (a), but BL occurs at the MT end for long 
MT. (d) a = 0.8, /3 = 0.9, u a = 0.001, ui d = 0.0008. There is no DW but BL occurs at the both end of MT for long MT. (a) (b) 
(c) are examples of l/(K + 1) < < 0.5, < /3 < 1/(A" + 1), 0.5 < j3 < 1 respectively, (d) is an example of 0.5 < a < 1 and 
0.5 < /3 < 1. 
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FIG. 2: Length of MT with high and low motor density, Nh and Nl, and their ratios Nh/N and Nl/N for three different 
cases in which DW occurs for large N limit, see Tab. |T] (a, b) a = 0.2, /? = 0.3, (c, d) a = 0.2, /3 = 0.7, (e, f ) a = 0.01, f3 = 0.1. 
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FIG. 3: The critical value of MT length N c decreases with attachment rates a,u a , and increases with detachment rates f3,ujd- 
Parameter values used in the calculations are (a) u) a = 0.001, uj d — 0.0002, (b) a = 0.2, /3 = 0.3. The meaning of critical value 
N c is that, DW occurs iff MT length N > N c . 
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FIG. 4: The velocity v dependence of motor traffic along MT. (a, b, c) are examples of steady state motor density for 
1/(K + 1) < /3 < 0.5, < P < 1/(K + 1), 0.5 < /3 < 1 respectively. The motor velocities used in (a) (c) are the same as given in 
(b). (d) are ratios of non-jammed MT length Nl/N, and (e) (f) are the jammed and non-jammed MT length (Nh and Nl) with 
the same parameters as given in (d). Other parameters used in the calculations of (a, b, c) are uj a — 0.001, aid = 0.0002, N = 150, 
and (a) a = 0.24,/? = 0.48, (b) a = = 0.1 (c) a = 0.25, = 0.9. The results indicate that N L ,N L /N increase with motor 
velocity v, and Nh,Nh /N decrease with v. 



